Proton Electrodynamics in Liquid Water 
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The dielectric spectrum of liquid water, 10'* — 10^^ Hz, is interpreted in terms of diffusion of 
charges, formed as a result of self-ionization of H2O molecules. This approach explains the Debye 
relaxation and the dc conductivity as two manifestations of this diffusion. The Debye relaxation 
is due to the charge diffusion with a fast recombination rate, 1/t2, while the dc conductivity is a 
manifestation of the diffusion with a much slower recombination rate, 1/ri . Applying a simple model 



based on Brownian-like diffusion, we find T2 ~ 10 s and ri 
the charge carriers, involved in each of the two processes, N2 - 



i 10 s, and the concentrations of 
5 X lO^*^ m-3 and iVi ~ 10^* m'^ 
Further, we relate and Ni to the total concentration of H30'''-0H~ pairs and to the pH index, 
respectively, and find the lifetime of a single water molecule, to — 10~^ s. Finally, we show that 
the high permittivity of water results mostly from flickering of separated charges, rather than from 
reorientations of intact molecular dipoles. 

PACS numbers: 77.22.-d, 66.10.Ed 



INTRODUCTION 

Electrical properties of water are of high importance in 
vital, environmental, and technological processes 
They have been under intensive investigation for 
many decades It has long been recognized that at 

room temperature, water is a good insulator with neg- 
ligible electronic conductivity and a dielectric constant 
e ~ 80. Potentiometric measurements reveal an appre- 
ciable proton conductivity, adc = 5-5 x 10~^ U~^m~^ at 
room temperature Jj]. This value is associated with the 
pH index, which is a key indicator of activity of protons in 
chemical reactions [sl, (9|. Normally, pH = 7; this water is 
regarded as neutral with the "free" -proton concentration 
of 6 X 10^^ m^''. It is believed that under normal condi- 
tions, a given H2O molecule will on average dissociate in 
roug hly 10"* seconds (11 hours) 

The origin of the high static dielectric constant of wa- 
ter is commonly explained by the orientational motion 
of the molecular dipoles, which is referenced as the De- 
bye relaxation The relaxation is particularly evi- 
dent in the frequency spectrum of dielectric permittivity, 
e*{uj) = e'{uj) + ie"{u!), as a strong anomaly near a char- 
acteristic frequency i^q — ujq/2tt; vq « 20 GHz at room 
temperature. 

Room-temperature dielectric spectra of water, re- 
plotted from the data of Refs. (j-lsl, are shown in Fig. 1. 
The s"{uj) spectrum is dominated by an absorption peak, 
accompanied by a step-like frequency dispersion in e'{u!). 
The step connects the high-frequency dielectric constant, 
£00 = 5, with the low- frequency (static) constant, e(0) = 
80. Let us note, that e'{uj) and e"(w) are connected by 
the Kramers-Kronig relations, thus the high value of the 



static permittivity, £(0), is due to a large integral inten- 
sity of the Debye absorption band in e"(w). 

The Debye relaxation in water, investigated experi- 
mentally and theoretically in many details 
is surprisingly well described by a simple relaxation for- 
mula: 



e'(w) = £00 + 
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Here td is the temperature-dependent relaxation time 
and Aej) — s{0) — e^o is the contribution of the dielectric 
relaxation to the static dielectric constant. 

For a long time, the Debye's idea about the orien- 
tational motion of water molecules has been exploited 
as the main microscopic mechanism responsible for the 
static permittivity. Basically since its introduction, it 
has been acknowledged, that the Debye model is over- 
simplified [H- Therefore, the model is being permanently 
modified. The present-day models involve the dynamics 
of protons and large molecular clusters Impor- 
tant is that in all these models, the geometry of the water 
molecule is a substantial input parameter. 

Here, we propose an interpretation of the dielectric 
spectrum of water alternative to the Debye's approach. 
We argue that considering exclusively the diffusive mo- 
tion of protons in water is sufficient for quantitative de- 
scription of its dielectric spectrum at frequencies lower 
than 10" Hz. 



REMARKS ON EXISTING EXPERIMENTAL 
DATA 

Due to the reasons, which will become apparent in 
the course of the article, our analysis is performed in 
terms of complex dynamical conductivity, rather than 
the dielectric function. The complex conductivity, a* = 
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FIG. 1: (Color online) Panoramic dielectric response of water 
at room temperature. Thin solid black lines are experimental 
data from Ref. . Dashed line is the frequency-independent 
(at < < 10^ Hz) conductivity according to Ref. 0, 
Thick solid magenta lines are calculations with Eq. [5] and the 
parameters from Table 1. Top panel: real part of the dielectric 
permittivity, £'{i^)- Bottom panel: dynamical conductivity, 
o{uj), left-hand scale, and dielectric losses, e" (uj), right-hand 
scale. 



Re{a*)+ilm{a*) , is related to the dielectric constant via: 
£*{u}) = Soo + icr* (uj) / {eqlo) , where £o is the free-space 
permittivity po| . Thus, the real part of conductivity, 
Re{a*) = (T, is merely the imaginary part of permittivity 
multiplied by frequency: <t(w) = weo£"(a;). 

The room-temperature cr(w) spectrum of water is pre- 
sented in the main frame of Fig. 1 (left-hand scale). In 
order to be specific, we use the data from Refs. 0, 01 for 
the dielectric constant, Aed, the relaxation time, tjj, and 
the dc conductivity, adc- The values of these parameters 
for 25 °C are listed in Table 1. 

The most important, for further consideration, fea- 
tures of the conductivity spectrum are the two well- 
distinguished plateaus, cti and (T2, situated at 10^ — 10^ 
and 10^° — 10^^ Hz, respectively. The plateaus are con- 
nected by a section, where a{io) oc lo" with s = 2. At fre- 
quencies around the lower-frequency plateau, the avail- 
able experimental data differ from each other. On the one 
hand, there are spectroscopic indications of a frequency- 
dependent conductivity in this range 

^ (th in black line 



TABLE I: Experimental data from Refs. [1, 0] used as input 
parameters in our model. All numbers are at room tempera- 
ture. 
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0-2 ~ 0-0(0; -> 00) 
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8.3 X 10"^^ s 


5.5 X 10"'^ n^^m"^ 


78 n^^m"^ 



in Fig. 1). On the other hand, in electrochemistry 
the conductivity below 10^ Hz is considered frequency- 
independent and equal to its dc value (dashed line 
in Fig. 1). Noteworthily, in the common conductome- 
try measurements, the working frequency is usually not 
even mentioned [7[. Hereafter, we take adc as the lower- 
frequency plateau value. Taking the data from Ref. Q 
instead, would only lead to minor quantitative correc- 
tions. 



ARGUMENTS FOR OUR MODEL 

Being expressed in terms of c(i^), the Debye relaxation 
form (Eq. [T]) looks like: 



(2) 



At high frequencies, auluj) is frequency-independent, 
corresponding to the second plateau of the experimen- 
tal spectrum: 



(72 « (T£)(t^ 00) = AeD£o/TD- 



(3) 



Surely, there is no physical difference between the rep- 
resentations in terms of £*(w) and in terms of (t*(w). 
However, we believe the conductivity representation 
gives a tip for a fresh look on the dielectric spectrum. 
Whereas the bell-shaped relaxation in s"{uj) is intu- 
itively connected with the orientational motion of the 
H2O molecules, the cr(ijj) curve of Fig. [T] hints to an 
alternative mechanism, namely, to the proton diffusion. 

In fact, the (j{uj) spectrum, given by Eq. [2l is consis- 
tent with acceleration-less motion of a charge q with mass 
m in a parabolic potential, 1^9(2;) = kx'^/2. The equation 
of motion for this charge is: 7717 x +kx = qE, where n is 
the spring constant and 7 is the relaxation rate. Then, 
the mobility is = q/lrwy), the diffusion coefficient is 
D = ksT^/q^ and finally the conductivity is: 



[q^N{kBT)/Di^]uj'' 

l+[(fci3r)2/K2^2]^2' 



(4) 



which coincides in spectrum shape with Eq. [51 

Noteworthily, the diffusion of particles, interacting 
with attractive centers, reveals such D{uj), that gives rise 
to the same dispersion in conductivity as in our Eq. |4] 

The cr(w) spectrum, consisting of two plateaus and a 
a (X uj^ section in-between of them, is typical for mate- 
rials with high ionic conductivity, the superionics [2^ . 
The conductivity spectra of these materials have been 
studied in great details. Their common feature is a sup- 
pression of the dynamical conductivity at low frequencies 
due to localization of diffusing particles in the minima of 
lattice potential and/or due to interactions between the 
particles (22| - [25l |. 
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We believe, that in regard to its proton conductivity, 
water gives all reasons to be compared with superionics. 
In accordance with the modern concept, protons, H+, 
and hydroxyl ions, 0H~, are permanently generated (due 
to the self-ionization of H2O molecules 13, U, 3 12 1) 
and recombined in the volume of water. 

Since free protons in water are not observed, they are 
considered to localize after their birth on neighboring 
neutral H2O molecules (on femtosecond time scale). The 
excess proton converts the H2O molecule into a charged 
complex HaO^ with a positive charge , and leaves a 
"hole", 0H~, with a negative (twin) charge q~ . Subse- 
quently, by a relay-race manner the and q~ charges 
wander diffusively over H2O molecules until they meet 
each other and recombine to produce a neutral H2O. 
The favorable-unfavorable molecular configurations for 
the proton exchange are stochastically formed by the 
thermal molecular motion 



"dc-conductivity" 
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The diffusion paths of the separated charges, from their 
birth to recombination, are sketched in Fig. 2. The twin 
q~^ and q~ charges are always in Coulomb field of each 
other. Therefore, in their majority they do not go far 
from the places of their birth (point 1 in Fig. 2). In- 
stead, they recombine with their own twin partner in a 
close vicinity of their birth places the [area with charac- 
teristic size £, point 2 in Fig. 2]. Sometimes, however, 
the "twins" fail to meet each other and recombine with 
"foreign" partners on a much larger distance L (points 3 
and 4). Obviously, the most probable foreign partner is 
a partner from the first configuration sphere of the ion- 
ized molecules. Thus, effectively, there are two recombi- 
nation processes for charges in water, faster and slower. 
We believe, that (similarly to superionics) the character- 
istic lifetimes of the two processes reveal themselves as 
characteristic knees in the conductivity spectrum. 



THE MODEL 

The model outlined above can be described by a set 
of equations. Let us designate the concentration, the av- 
erage lifetime, and the mean free path of the charges, 
involved in the slow and in the fast processes, as A'^i, ti, 
L, and iV2, T2, i, respectively. A'^i and A'2, can also be in- 
terpreted as the concentrations of the HaO^-OH" pairs 
participating in the slow and fast processes. Similarly, 
Ti and T2 are the average life times of the HaO^-OH" 
pairs, and L and £ are their characteristic sizes. 

The basic assumption we made is that the lower and 
upper conductivity plateaus in Fig. 1, ai and a2, corre- 
spond to diffusion of long- and short-living charges, the 
latter (the fast) process giving rise to what is commonly 
refereed as the Debye relaxation. Hence, the complete de- 
scription of the a{uj) spectrum can be written as a sum 
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FIG. 2: (Color online) Schematic diagram of the two-scale 
proton conductivity in water. The black dots are protons; 
the shadowed small circles are oxygen atoms. The dashed 
lines show diffusion paths of the charges g"*" and q~ , the single 
dashes representing single proton hops. The large gray circles 
are coordination spheres of ionized molecules. For further 
explanations, see text. 



of two terms: 



a (w) = CTi 



1 + t^Vf 
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(5) 



The conductivity at each, the lower and the up- 
per, plateau is connected to the diffusion coefficient of 
charges, D, by the Nernst-Einstein relation: 



(Ti 2CNiD; 0-2 = 2CN2D, 



(6) 



where C ~ q^/ksT and the coefficient 2 takes into ac- 
count the occurrence of positive and negative charges. 

Using the Einstein-Smoluchowski formula, we can also 
connect L, ri, £, and T2 to the diffusion coefficient: 



D = 



- ill 
6t2 6ti 



(7) 



Because both, slow and fast, processes span over all 
volume of water, for the unit volume one can write: 



(8) 



1 = ^L^N, = ^£^N2. 
3 3 

Equations [S]- [5] constitute an equation set, the analyt 
ical solution of which gives the following result: 

2 



N2^3i4nf{^y;£ 



47rT2CT2 6T2 v47r<72 



(9) 



Ni^ N2X—;L^£x ( — 



1/3 / \ 2/3 



; Tl = T2 X 



(10) 

The numerical values of these parameters, calculated 
using the experimental data of Table 1, are presented in 
Table 2. With these parameter values, Eq. [5] compre- 
hensively describes both, the Debye relaxation and the 
dc conductivity (thick magenta line in Fig. [T|). 
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TABLE II: Room-temperature numerical values of the param- 
eters obtained from our model (Eqs. [91- I10p . 



iVi 


N2 


D 


3.8 X 10^^ m-=* 


5.4 X 10^'^ m"^ 


1.2 X 10"* mVs 


L 


£ 


Tl 


0.18 fira 


0.76 nm 


5 X lO"'^ s 



The found parameters characterize a Brownian motion 
of thermally activated charges and q~ over the "sea" 
of neutral molecules, i.e. the drift currents. The mutual 
thermal motion of neutral H2O molecules is not relevant 
for our consideration. 



CONSEQUENCES OF THE MODEL AND 
CONCLUSIONS 

1. Concentration of separated charges and life time 
of H2O molecules. As one can see from Table 2, the 
concentration of short-living H30"'"-0H~ pairs is huge, 
N2 0^ 5 X 10^^ m-3. Because N2 > A^i, N2 can be 
taken as the total concentration of the HaO^-OH" pairs. 
Thus, roughly 1% of all H2O molecules are ionized, the 
concentration of H2O molecules being TVg = 3 x 10^^ 
m~'^. The found charge concentration is by several or- 
ders of magnitude larger than the commonly accepted 
value for neutral water, 10~^ mole/liter = 6 x 10^^ m~'^. 
Figuratively, water constantly boils with ^-sized HsO"^- 
0H~ pairs, which have a life time of T2 ~ 10 ps. This 
result agrees with what was reasoned in Ref. P^] based 
on molecular-dynamics stimulations. 

Our value for the size of HsO^-OH^ pairs, £ — 0.76 
nm, correlates with the recent x-ray scattering results 
of Ref. (26l |. where density fluctuations in water were 
reportedly found on a comparable scale. 

Because A^2 >> ^1, the life time of a neural H2O 
molecule in thermodynamic equilibrium, tq, can be esti- 
mated from the following equation: 



N2 

T2 



No 

To ' 



(11) 



From here, we obtain tq ~ 1 ns, which is 14 orders of 
magnitude smaller than the "standard" 10*^ s. 

2. Permittivity. The Debye relaxation time changes 
its meaning - it is now the average life time of the short- 
living separated charges (the HaO^-OH" pairs). For the 
permittivity spectrum, by combining Eqs. |31 [H and HI 
we obtain: 



Ae 



D 



(g X e)^ N2 
3kT eo 



(12) 



i.e., the pairs of the separated charges can be considered 
as dipoles with the average dipole moment p = q x i. 



The polarizability (per unit volume) of these dipoles is 
A^2 X p^/ikT. Thus, Eq. [T2l shows that the diffusion of 
the charges to the average distance t is the reason for 
the step in the dielectric function, Aed, at frequencies 
around l/27rT = 10^" Hz. The experimental value of this 
step, A££) « 75, is automatically fulfilled in our model 
(through A'2 and t). This means that the dominant con- 
tribution (more than 90%) to the static dielectric con- 
stant (e(0) « 80) is provided by the restricted-distance 
currents of separated charges, rather than by orienta- 
tional relaxation of intact H2O dipoles. 

Obviously, the orientational motion of the intact 
dipoles should also reveal itself in the dielectric spectrum. 
We believe, a good candidate for this is the Debye-like 
bands, found in the measurements at frequencies higher 
than 10" Hz 
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28l | . These higher-frequency relaxation 
processes are also seen as the upturn of the experimental 
conductivity at the highest frequencies in Fig. 1. 

The slow recombination process with the characteristic 
time Tl also provides a contribution to the static dielec- 
tric constant. According to Eq. [3J this contribution is 
however very small, Aei < 1. 

3. Diffusion coefficient. In our consideration, the dif- 
fusion coefhcient, D ~ 1.2 x 10^^ m^/s is related to the 
relay-race diffusion of charges ((7+ and g~), not to diffu- 
sion of a tagged proton. In literature, however, this value 
(more accurately, 9.3 x 10^^ m^/s) is commonly accepted 
as the diffusion coefficient of protons [§| and considered to 
be "anomalously" high. According to our findings, the 
real Brownian diffusion coefficient of a tagged proton, 
D"+, is 100 times smaller: D"+ = {a^/&To) ~ lO'^" 
m^/s, where a = 2^ A is the distance between the H2O 



molecule centers [17 1. 



It is worth noting here, that if one substitutes the pro- 
ton diffusion coefficient, D^^ , and the total proton con- 
centration, 2A^o, into the Nernst- Einstein equation (Eq. 
[6]), one gets the upper limit for the proton conductivity 
in water: cr = ACNoD"+ = 74 fi-^m-^ that is pre- 
cisely the value of the higher-frequency plateau in the 
o{ui) spectrum, (T2. 

4. dc conductivity and relevance to pH index. The 
found value for Ni (3.8 x 10^^ m~'^) is practically equal 
to the commonly accepted concentration of "free" pro- 
tons (6 X 10^^ m~^), which provides pH = 7 in neutral 
water. This result shows that the common dc conducto- 
metric methods detect only those "survived" protons (of 
concentration Ni ) , which are involved in the slow recom- 
bination process, while the short-living protons (concen- 
tration A^2) are not detectable in these measurements. 
The occurrence of the short-living protons in water re- 
quires introducing a new, "fast", dissociation constant, 
Kw2^ in addition to the common ("slow") Kyyi, related 
to the pH index. Since both, long- and short-living, pro- 
tons are chemically active, the conventional conception 
of the pH index could probably be revised in such a way 
that the fast Kw2 is also taken into account. 
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Summarizing, we have found that the dielectric spec- 
trum of hquid water at frequencies below 10^^ Hz can be 
entirely understood in terms of proton diffusion, resulting 
from self-dissociation of H2O molecules. No long-living 
geometric structures, created by the water molecules, are 
needed to be included in this consideration. We believe 
that many other physical properties of water are also de- 
termined mostly by the ability of H2 O molecules to dis- 
sociate rather than to form any sorts of geometric struc- 
tures. 

We are grateful to S. D. Zakharov and G. M. 
Zhidomirov for useful discussions. 
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